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We discuss resummation strategies for free energy in quantum field theories at nonzero temper- 
atures T. We point out that resummations should be performed for the short- and long-distance 
parts separately in order to avoid spurious interference effects and double-counting. We then discuss 
and perform Pade resummations of these two parts for QCD at high T. The resummed results are 
almost invariant under variation of the renormalization and factorization scales. We perform the 
OA ■ analysis also in the case of the massless scalar <j> 4 theory. 
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I. INTRODUCTION 



During the last decade the physics of many-particle systems at relativistic energies has attracted increasing interest. 
On the one hand, this was due to the hope that high-energy ion collisions at present day colliders would reach the 
kinematical regime of the quark-gluon-plasma phase of QCD and would thus give new insight into the phase structure 
j of QCD and into the confinement phenomena. On the other hand, a better understanding of such systems can be 
Q\ • applied to ultrarelativistic astrophysical and/or cosmological configurations and should therefore provide a conceptual 
tool for describing the physics of the early universe, in particular phenomena such as baryon number generation, big- 
bang nucleosynthesis, etc. Although it is not yet clear to what extent the various mentioned configurations can be 
considered as being in thermal equilibrium, the first attempts at a systematic analysis concentrated on equilibrium 
statistical mechanics. 

Since we are convinced that the interaction of particles at energies up to 1 TeV is described by the standard 
model of strong and electroweak interactions, the main task consists in formulating a consistent thermal quantum 
field theory of QCD and QED. Because QCD is asymptotically free, the corresponding running coupling parameter 
O ,' becomes small at a sufficiently high temperature T (a s (T) <C l) 1 and a perturbative treatment seems possible. This 
expectation is supported by the observation that lattice gauge theory calculations jl], |^, [| indicate that the physics 
1 of the (quark-)gluon plasma at T larger than four times the critical temperature T c is to a very good approximation 
that of an ideal gas. In particular, the free energy F of such a system is off the ideal gas value by less than 20%. Even 
more, if masses with an appropriate temperature dependence are accepted, the lattice results [|4| agree almost exactly 
with the ideal gas ones for all values of T > T c . This would suggest that perturbative treatments remain candidates 
jj^ ■ for a description of the relevant phenomena. 

Unfortunately, a straightforward perturbative treatment has led to results which are far from satisfactory. The 
problems emerging can be seen best in the expressions for the free energy F as calculated in various theoretical 
models. The most exhaustive calculations, up to ^g 5 , have been performed for the massless </> 4 theory [| and for 
QCD (without and with quarks) jj], ||. These results have the following unpleasant features: 

1. There is a very bad convergence behavior of the series, 2 some of the higher-order terms in the perturbation 
series are larger than the leading term unless the coupling parameter is very small - in QCD, the coupling a s (T) 
would have to be at T > 1 TeV. If successive orders are included in the truncated perturbation series (TPS), 
the sum takes on alternating large values (cf. Fig. 1 of Ref. S). 

2. At every fixed order the TPS's show a strong dependence of the (arbitrary) renormalization scale (RScl) /x. This 
results in an additional severe uncertainty of the evaluation, this time due to the ambiguity of the choice of the 
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T being considered as a measure of the average energy of the system's constituents. 
2 Formally, the series is believed to be an asymptotically divergent series. 
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RScl fi. In an asymptotically free theory one usually is on the safe side if one chooses \i close to the (lowest) 
physically relevant momentum scale typical of the contributing configurations, since otherwise the expansion 
coefficients would blow up. In the case of a system at temperature T one would be led to choose \i ~ T, or 
more specifically /i — 2t:T, i.e., the energy of the first nonvanishing Matsubara mode. In fact, this is the choice 
actually used in the literature. However, we should keep in mind that in the case of a plasma consisting of 
massless components there are important physical effects connected with much lower energies (screening of the 
chromoelectric or chromomagnetic field). Therefore, the choice fi ~ 2ttT seems not to be well founded, a priori. 

These features are related to the rather large genuine collective effects. And, in fact, it has been noted recently 
that these unexpectedly large corrections are rooted in the terms ~g 3 ,g 5 within the expansions. Such terms render 
the expression nonanalytic in a = g 2 /(47r 2 ) and they do not emerge in the ordinary (T = 0) field theory. They are the 



result of a certain (partial) resummation which is necessary to get rid of these infrared divergences [O, Q . They are 
closely related with collective effects such as screening or Landau damping. 

This observation has constituted the basis for several attempts to remedy the situation: since a specific method 
of resummation seems to be the root of the problems, one is led to try different resummation procedures. The main 
attempt goes under the name of screened perturbation theory (SPT) |h], |l^| and has been made more systematic by 
what is called optimized perturbation theory Jl3| . The main idea is to add a local mass term to the free part of the 
Lagrangian and subtract it from the interaction part. The former is treated nonperturbatively, constituting a genuine 
(screening) mass; the latter perturbatively. Physically this means that one expands about a (ideal) gas of massive 
quasiparticles (dressed gluons, for example) rather than about massless particles. The technical consequences are 
striking. Since the objects one starts with are massive, there is no infrared problem and no need for resummation, 
and the resulting expressions have better convergence behavior than the original perturbation series. However, for 
obtaining numerical results, one has to fix the chosen mass at the final stage. The authors of Ref. p% , e.g., used (an 
appropriate approximation of) the gap equation for fixing this mass. 

All this is relatively straightforward in a theory such as cf> 4 where nothing forbids a genuine mass term. Within a 
gauge theory like QCD the whole procedure becomes much more cumbersome, since the addition of a genuine (local) 
mass term would spoil the gauge symmetry from the outset. A technical way out of this problem is hard-thermal 
loop (HTL) perturbation theory |J. This is a SPT generalization which respects gauge invariance. The procedure 
again rests on a (partial) resummation: Those higher-order loop corrections are included that are of leading order in 
g s for amplitudes involving soft external momenta p ~ g s T. Unfortunately, the resulting HTL correction terms are 
nonlocal. This leads to complicated UV divergences, and only some of them are canceled by physical mechanisms 
(quasiparticle formation and Landau damping). The rest have to be tamed by artificial counterterms proportional to 
the quasiparticle mass and they generate an additional renormalization scale dependence. Despite this complicated 
and not uniquely specified procedure, the resulting series has again a better convergence behavior than the original 
perturbation series. On the other hand, when the so called <&-derivable approximation scheme jl4|, is applied, 
some of the problems of the HTL perturbation theory are avoided. In this scheme, the HTL contributions to the free 
energy were resummed in Ref. |D| in 4 theory, QED, and QCD, but the UV divergences were shown to cancel. A 
related approach which uses the framework of the scalar 0(7V)-symmetric model in the large- N limit was developed 
in Ref. |lq| . Another approach, overcoming the infrared problems appearing in the conventional perturbation theory 
and using dressed propagators, was developed in Ref. J17[. 



On the other hand, the strong renormalization scale(RScl) dependence is not or not sufficiently reduced by the 
SPT or HTL approach. Further, there remains a principal question: Why are the SPT results so much better than 
those of the naive perturbation approach? Although the need for resummation is avoided in the first step in SPT 
due to the mass m g of the quasiparticle, the correct choice of the specific formula for m g again needs a resummation 
implicitly contained in the gap equation. As long as these points are not clarified, the SPT treatment cannot be 
considered fully satisfying. 

A completely different way of remedying the weak points of the naive perturbation theory consists in replacing 
the truncated perturbation series (TPS) by appropriate Pade approximants (PA's). This approach is motivated 
by at least two features of PA's: First, PA's at increasing order in general show much better convergence than 
the TPS's from which they are obtained. It has been shown that even when the TPS's are divergent (asymptotic 
series), the corresponding PA's may converge and do it under rather general conditions Jl~8| | . Secondly, PA's reduce 
considerably the RScl dependence of the TPS. In fact, it is known that the diagonal PA's, constructed from a 
TPS of an RScl-independent quantity in powers of a(/i) = g 2 (/i)/(47r 2 ), are RScl-independent in the limit of the 
one-loop running of the coupling parameter a(/x) (large /3o-lmrit) InJ]. Further, related approximants have been 
developed which are exactly RScl-independent pfj| and even renormalization scheme independent pl| . Since the 
full perturbation series corresponds to a physical (in principle measurable) quantity, this quantity is exactly RScl- 
independent. Therefore, one is led to conjecture that PA's and related approximants are nearer to the true value 
than the original TPS's. The physical reason for the (approximate) RScl independence of such approximants is that 
they include a certain resummation [[L9| , thus containing infinitely many terms whose absence was responsible for the 
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spurious RScl-dependence of the TPS. 

The first applications of PA's to thermal perturbation theory were made by Kastening pi| and by Hatsuda |p3| . 
These authors started from the the available TPS for the free energy F (in (j) 4 and QCD), which is a TPS in powers of 
a 1 / 2 = g/(2w), up to ^a 5 ^ 2 , and replaced it by various PA's based on it. They demonstrated improved RScl stability 
of these PA's. Although the relatively large number of powers in g of the underlying TPS suggests a good convergence 
quality of the resulting PA's, the results have to be treated with caution. First, there is no formal reason to expect 
that certain PA's, specifically the diagonal ones, should be stable under variation of the RScl. In this respect, we note 
that it is the one-loop (large- /?o) running of a = g 2 /(47r 2 ), 

1 + a(/z )po M/XjlM)) 

which is responsible for the (large-/?o) RScl independence in diagonal PA's of underlying TPS's in powers of a. This 
is so because relation ([j]) represents a homographic transformation a i— > a/(l + Ka) |TSt| . However, relation ([j]) yields 
a structurally different one- loop running of g = 2ira 1 ^ 2 , namely 

This relation is not homographic, and consequently the PA's in such a case need not be (large-po) RScl-invariant. 

There is a second, more significant point of criticism of the aforementioned PA approach. It can be best explained 
in terms of diagrams. As we have mentioned, the PA's represent a resummation (analytic continuation) of the 
infinite sum of a certain class of diagrams. On the other hand, the TPS expressions for F(T) also include a selective 
resummation (of ring and super-ring diagrams). Therefore, a naive application of PA's to the original TPS constitutes 
a mixing of (at least) two inequivalent resummation effects and could easily lead to partial double-counting or other 
(interference-like) inconsistencies. In the context of QED and QCD at T = 0, somewhat similar aspects have been 
pointed out and accounted for in Refs. j24|. 

Within the present paper we want to stick to PA's because of their unique advantages: reduced RScl-dependence, 
better convergence properties, and (quasi) analytic continuation aspects. However, we are going to develop a procedure 
that is free of the aforementioned weaknesses. We do this by separating all the terms in the full available TPS into 
groups which represent TPS's of separate physical, i.e., RScl-independent, quantities. To each of the TPS's we 
apply PA's separately. The resulting expressions are not only (approximately) RScl-independent, but are supposed 
to represent better approximations to the true values since double-counting and interference-like inconsistencies are 
excluded. 

In Sec. II we present the main elements of our method. In Sec. Ill we apply this method to QCD, and in Sec. IV 
to the massless 4> 4 theory. Section V summarizes the results and presents conclusions. 

II. SEPARATION OF LONG AND SHORT DISTANCE REGIMES 

We start with the perturbatively calculated expressions for the free energy density F both for massless <fi and for 
QCD. They have been calculated recently ||, ||, 0, || up to ~g 5 , i.e., three loops plus all the ring diagrams summed 
up. The generic structure of the resulting expressions is 

F(T) = F idca i [C + C 2 a + C 3 a 3/2 + C 4 (n)a 2 + C 5 (fj)a 5 / 2 + 0(a 3 In a)] . (3) 

Here, a is short-hand for a(/i) = g 2 (/i)/(47r 2 ), i.e., the running (RScl-dependent) coupling parameter. The renormal- 
ization scheme (RSch) is assumed to be fixed, say MS. The RScl-running is described by the perturbatively specified 
renormalization group equation (RGE) 

|^4 = - / 3oa 2 -/3 1 a 3 -/3 3 a 4 -... . (4) 

It is interesting to note that in the series (||) the coefficients Co,C2,C3 do not depend on the RScl n, but C4 and 
C5 do, the /i-dependence showing up as an additive contribution oc ln[/i/(47rT)] to the corresponding coefficient. In 
QCD, C4 includes a term ocln[a(/i)] [~ 1/ lnln(/i)]; in </> 4 theory, C5 includes a term ocln[a(/z)]. It is crucial to keep in 
mind the origin of the various terms: The terms ~ a 3 / 2 and ~ a 5 / 2 come exclusively from the resummation of the ring 
diagrams that is necessary to avoid infrared divergences g, f7|. This resummation procedure also yields part of the 
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coefficient C4 of the ~a 2 term. In particular, the contributions ocln(a) in C4 and C5 are generated by this procedure 
as the remnant of the emerging dependence on the screening mass. The other terms (Co, C2, and the remaining part 
of C4) stem from ordinary perturbation theory in powers of a and do not contain resummation effects. We further 
note that the term ^a 3 (~g e ) cannot be obtained by perturbative methods (which include ring summation) because 
of the severe infrared divergences appearing at that order p5| . 

Our aim is to apply Pade approximants (PA's) in a consistent manner to a TPS of the type (H). As argued 
in the Introduction, in order to avoid an uncontrollable mixing and superposition of different resummations, the 
separation of the pure perturbative from the ring resummation-generated terms should be performed. In addition, 
however, we want to take advantage of the approximate /^-independence of the PA's when they are applied to TPS's 
(in a) of /i-independent quantities. Therefore, the right-hand side of Eq. (||) should be split so that the resulting 
parts (infinite power series) represent quantities that are /i-independent ("physical") separately. This suggests that a 
physical principle of separation should be involved, i.e., one which is connected with measurable effects. 

In finding such a separation principle one is guided by the decomposition of all thermodynamic quantities in 
(Fourier) modes. Within the imaginary time formalism, the free energy F at every given order is expressed as a sum 
over the Matsubara frequencies u> n 

uj n = 2nTn for bosons , (5) 
= 7rT(2n + 1) for fermions , (6) 

where n = 0, 1, . . .. Since the thermodynamic quantities and correlation functions can be derived from F, they also 
show up as sums over modes. For a given correlation function the contribution from the (exchange of the) Fourier mode 
with frequency ui n falls off at large spatial distances as exp(— uj n R) (at least if T is larger than all contributing masses). 
Therefore, the only mode which does not produce an exponentially vanishing contribution to the long-range correlation 
functions is the bosonic zero mode u>q = [Eq. (^)] . Consequently, at sufficiently high T, the static correlators of the 
contributing fields at large distances R 3> 1/T are exclusively determined by the zero mode. Since this long-distance 
behavior of correlators is, at least in principle, observable ("physical"), the procedure of separating the bosonic zero- 
mode (long-distance) contributions from all the other (short-distance) ones rests on physical grounds. We therefore 
expect that both contributions are separately /i-independent since both have a physical meaning. Further, we know 
that all resummation effects in the series (|^) contribute exclusively to the long-distance part. This is so because these 
resummation effects of the ring diagrams exclusively stem from the zero- mode contributions, which represent the 
strongest infrared divergences at the single diagram level. The long-distance part, i.e., the rcsummed ring diagrams, 
shows up as a power series in powers of g oc a 1 / 2 , starting with the a 3//2 term. The short-distance part, on the other 
hand, has the ordinary perturbation character - a perturbation series in integer powers of a = g 2 / (An 2 ) . 

The discussed decomposition of F represents the basis for our improvement of the underlying TPS results: We 
apply the appropriate PA's to the separate parts, so that (at least some of) the approximants are approximately 
RScl- independent, and presumably better converging. In the case of the short-range contributions these are the 
diagonal PA's [n/n](a). At the available order of the underlying TPS (|J), the short-range part is of the form 
Fg(T) = Fideai(Co + C%a + C4C1 2 ), and the only possible diagonal PA is [l/l](a). 3 In the case of the long-range 
contributions, one tries to see which PA [n/m](g) (n + m = 2) is approximately RScl-independent. 

Before continuing, we mention an alternative way to interprete the described decomposition on physical grounds: 
The long-distance contribution reflects the observable phenomena of screening of interactions and formation of quasi- 
particles (collective modes). To understand this, we recall the following facts. It is with the resummation of the ring 
diagrams that the zero modes acquire a (screening) mass m s ~ gT (m s can be defined as the solution of the gap equa- 
tion 0^). And vice- versa: It is the screening phenomenon which explains the striking fact that, while all individual 



ring diagrams are infrared divergent, their total sum is convergent. From this one concludes that, due to the observable 
nature of the screening phenomenon, the long-range part Fl, which encompasses the resummation-generated effects, 
is of physical nature and should therefore be RScl-independent. This can best be seen by determining the screening 
mass via the gap equation |l^, [T^] which is equivalent to summing up the ring diagrams. 

This interpretation also shows that in theories with more degrees of freedom (like QCD) there might be several 
stages of screening, e.g., screening of chromoelectric and chromomagnetic gluons, respectively. Consequently, the free 
energy F can then be decomposed into more than two parts: the short-distance part, and one part for each kind of 
screening. 

Having described the physical idea behind our approach, we have not yet addressed the technical problem of actual 
decomposition. In principle, the answer is simple: in order to single out the long-range (zero-mode) contribution, one 



3 The Pade approximant [n/m] (x) for a quantity S(x) is a ratio of two polynomials of ra'th and m'th degree in x: [n/m] (x) = P n (x)/ P m (x), 
such that its expansion in powers of x reproduces the terms up to ~x m + n of the expansion of S(x): S(x) = ro+r\x+. . . r n + m x n+m +. . .. 
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has to integrate out all higher modes (bosonic, and all fermionic modes). The resulting expression, when subtracted 
from the full expression, should yield the short-distance contribution, i.e., that of all the nonzero modes. In practice, 
integrating out explicitly all higher modes is a cumbersome task. Fortunately, there exists an alternative method, 
the method of effective field theories, which was developed for thermal perturbation theory by Braaten and Nieto 
H H |2(|. The method can be briefly described as follows. At small distances (R < 1/T) the behavior of the system is 
determined by ordinary perturbative QCD. On the other hand, the long-distance behavior is dominated by the zero 
modes and can, consequently, be described by an effective field theory which, at large T, is a bosonic field theory 
in three dimensions ("dimensional reduction" p5| , ^7|). The effective bosonic field is approximately identified with 
the zero modes of the original fields. For the construction of the effective field theory one does not need to specify 
this effective (static) boson field in terms of the original fields exactly. One simply writes down the most general 
three-dimensional Lagrangian for the effective fields that respects the symmetries of the original theory. In general, 
this effective Lagrangian contains infinitely many terms - operator expressions of arbitrary high dimensions - and is 
thus nonrenormalizable. Therefore, an ultraviolet cutoff A is needed. The corresponding effective coupling parameters 
9E,i (i=lj 2, 3; E stands for effective) are then determined by a matching procedure: One computes sufficiently many 
static correlation functions, both in the original and in the effective theory, with as yet unspecified gE,i, and demands 
that the results agree at larger distances R > 1/T. The resulting matching relations allow one to express the effective 
gE,i's in terms of the original coupling parameter g, T, and the cutoff A: 



The cutoff A (~ T^2ttT), also called the factorization scale, is roughly the momentum below which the effective theory 
should take over. Relations (0) can be understood as perturbation series in the coupling parameters g. Fortunately, 
the g_E,i's corresponding to interaction operators with higher dimensions are of higher order in g. Therefore, for 
calculation at a given order only a very restricted number of gE.i's has to be taken into account. 

The effective theory approach has been applied by Braaten and Nieto to 4 theory || and to QCD ||. 

In the case of <7 2 4 /4! there is only one cutoff A separating the long- and short-distance regimes ||. Consequently, 
the free energy density F consists of two parts 



where the long-range part FL is determined by the effective theory and contains the effective coupling parameters 
gE,i- If these gE,i's are expanded in powers of g, then Eq. (ph yields the formula for F that was originally obtained 
perturbatively by including the ring-diagram resummation [t]] . It is a check of consistency of the effective field 
approach to show that the total expression for F is independent of the factorization scale A. 

In QCD one can separate from the short-distance regime two long-distance regimes (chromoelectric and chromo- 
magnetic) || corresponding to two types of screening. Therefore one can apply two effective theories. One is called 
the electrostatic QCD (EQCD), and contains electrostatic and magnetostatic effective gauge fields. It describes the 
physics of the quark-gluon system at distances r ~ I /toe, where wie ~ gT is the mass scale of the chromomag- 
netic screening. The other is called magnetostatic QCD (MQCD), and contains only magnetostatic fields. It acts 

at distances r ~ 1/tom, where tom ~ g 2 T is the magnetic screening mass, which, however, cannot be calculated 
perturbatively because of severe infrared divergences at the order g 6 [^5|, |2q| . As a consequence, the free energy 
density F of hot QCD consists of three contributions: 

1. the short-distance [r ~ I/(-7rT)] contribution Fe of nonzero modes, i.e., the modes with frequencies equal to and 
higher than the first Matsubara frequency lo\ = 2'kT\ 

2. the long-distance contribution Fm, with r ~ 1/toe ~ l/(gT), i.e., the collective modes described by EQCD 
effective theory; and 

3. the "rest" contribution Fr from even larger scales r ~ 1/mM ~ l/(# 2 ^) °f MQCD effective theory. 
Each of these three contributions is expected to be physical, i.e., fi- independent. 



9E,i = gE,i{g, T, A) . 



(7) 



F = F s + F L , 



(8) 



III. THE CASE OF QCD 



A. Formulas for the approach 



We will first apply our approach to QCD since this is physically and experimentally the most interesting case. The 
free energy density F, for an arbitrary number rif of quark generations, has been calculated in thermal perturbation 
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theory up to terms ~a 5 / 2 [Q, 0), within the MS renormalization scheme. As mentioned before, the free energy density 
can be decomposed into three physically distinct parts |p| 



F — Fe + Fm + F] 



(9) 



Up to order a 5 / 2 , only Fe and Fm contribute, and they are, in principle, perturbatively computable. These are the 
contributions from the energy intervals (Ae, oo) and (Am, Ae), respectively, where the factorization scales Ae and Am 
delimit the ultraviolet (UV) and infrared (IR) bounds of EQCD and satisfy the inequalities 



mM(~g 2 s T) ~ A M ~ m E {~g s T) ~ A E ~ wi(=27rT). 



(10) 



Apart from fitting into this hierarchy, the factorization scales Ae and Am are in principle arbitrary. The sum of the 
three contributions is independent of these scales, and Fe + Fm is independent of Ae- Once Ae and Am are fixed, 
each of the three contributions is a (quasi)observable, i.e., a quantity that is independent of the renormalization scale 
and scheme. We will show that the specific RScl-dependence of the available TPS's for Fe, tue, and Fm is in fact 
consistent with RScl-independence of the full quantities Fe, the, and Fm- 
Specifically, we have for Fe the following TPS flcj: 



Fe 



87T 2 

45 



21 \ 15 , 

1+ 32»' -T" 



12 



/ ^e(Ae) 



where the first term represents the contribution of the free (ideal) quark-gluon gas 

21 



8tt 2 4 

r ideal a r- 

45 



1- 



32 



and the "canonical" QCD-part Fe(A e ) is 
-Fe(Ae) = o,(fi 2 ) 1 1 + 



Q(n f ; A E ) + L(n f ) In (j^J 



a{^)}+0{a A ) , 



with 



Q(«/;A E ) 
L(n f ) 



= -18 



( 1 + "//6) ln f A 2 



(l+5n//12) 
11 f } _2n f 



V 4tt 2 T 2 



4 (-244.898- 17.2419 ^+0.415029^) 
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(l+5n//12) 



33 



(11) 
(12) 

(13) 

(14) 
(15) 



Here, rif is the number of active quark flavors; a(fj, ) = a s (fi ;MS)/n = g 2 (/i , MS)/(47r ). The RScl /i is usually 
chosen as ^irT. We note that L(n/), the coefficient of the (In //^-dependent part at order a 2 , is exactly equal to f3o, 
the one-loop coefficient of the QCD beta function of the RGE (Q). It is then exactly this fact which guarantees that 
the derivative dFE/dln fj, 2 has the terms of ~ a 2 canceled due to the RGE. Therefore, this variation is ~a 3 , which 
is of the order of the first unknown term in the TPS of Fe. This shows that the TPS behavior of Fe (and Fe) is 
compatible with the supposition of the /i-independence of Fe- 

Further, the available TPS for the electric Debye-screening mass is S 



4tt 2 T 2 (l+n//6) 



a(n 2 



P(n/)+ /?o In 



4tt 2 T 2 



a(fi 2 



0(a J 



where 



P(n f ) = (0.612377- 0.488058 n f - 0.0427979 n 2 ) /(l + n//6) . 
Again, it is possible to see that Owle/ d\n [i 2 ~ a 3 , compatible with the /i-independence of this quantity. 



(16) 



(17) 



4 We note that there are two misprints in a formula for F E (/ E = F E /T) in Refs. (| - in Eq. (7) (PRL) and Eq. (54) (PRD) - in the 
sign of the coefficient at ln[A E / (2irT)] and in the sign of the term 17.24 nj. 
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Let us now consider the part Fvi, which includes all long-distance contributions, i.e., all ring-diagram resummation 
effects. The available TPS for Fm was calculated within EQCD (with effective mass toe and effective coupling <je) in 
Ref. M. Since <? E (/i) = g 2 (p)T [l+0(g 2 )], one can express the result immediately in powers of g s . This gives 



*M = ^-Tml F M (A E ; A M ) , 

07T 



(18) 



with 



A 2 

I 

2 



in 



F M (A E ;A M ) = -1 + (0.255838 + In I — (l + nf /6)-^ 2 g s (ii 2 ) 



2tt 



+27.5569^(1 +n / /6)- 1 <? s V) + 0{g 3 ) . 



(19) 



The AM-dependence would first show up at terms of higher order ~ g 3 . The A E -dependence cancels in the sum Fe+Fvi 
to the orders available in Eqs. (0)-(f|) and (|l|)-(|l|). Equation (H) is a series m g s (/j 2 ) = 2ir^a{i± 2 ). To the order 
available, it is automatically compatible with /^-independence of Fm since dg s (ii 2 ) / '9 In /j 2 r 
When the screening mass toe ( |l6|) is expanded in powers of g s {n) and inserted in Eqs. (|H 



g%hy the RGE 



19|) , a power expansion 



of Fm in g s (fJ-) can be obtained, starting with g 3 



777. 



3tt 



-k 3 T 4 



47T 2 fc 2 



fP(n/)+A)ln( 



E) 

0.255838 

2 



^ln(. 9s ( M 2 )fc / )+^ln(A E /T) 



A* 



4tt 2 T 2 



fc 2 + 27.5569 



<?,V) 
g 5 s (p 2 )+0(g 6 s h-ig s )} , 



(20) 



where we used the notation kf = (l + rij/6) 1 / 2 and Eq. (|T7|). Adding it to the expansion (pd|)-(|l3|) for Fe, a TPS 
for Fe+m in powers of g s (u) up to ~ g% is obtained as given in Ref. [[7| (third entry) and in Ref . g (second entry), 
i.e., a TPS of the form (||). Interestingly, due to the use of the expansion of toe, the coefficient Ci(p) becomes 
<7 s -dependent [dependent on lng 2 (/i)], which represents an additional, although possibly only formal, obstacle for the 
direct application of PA's to such a TPS. 5 We avoid this problem by using in the free energy the Pade-resummed 
squared screening mass m E . 

It is well known that the TPS @ for F M in QCD has very bad divergent behavior, and this is the case to a 
somewhat lesser extent for the TPS of toe. Thus, their direct evaluations do not yield useful predictions. Of course, 
this bad divergent behavior is transported to Fm and to F = Fe+Fm. To remedy this, the authors of Refs. |2^ |23] | 
presented evaluations of F via various Pade approximants (PA's), which were based on the TPS of the expansion of 
F in powers of g s , i.e., the sum of Eqs. (|Tl|)-([l3|) and (p0|), and no separation was performed. Although their results 
showed significantly reduced /7-dependence of the Pade resummed values of Fe+m in comparison to the /x-dependence 
of the TPS of Fe+m, we believe that this approach is not well motivated. This is so because it probably leads to 
partial double-counting or other interference-like inconsistencies, as argued in the previous sections. We will illustrate 
this argument with the following simple example. Suppose that we want to resum, by a PA, the sum S = (S1+S2) of 
two observables Si , S2 . Suppose that we have these two observables available as power series of a up to next-to- leading 
order (NLO): 



Sj = a(l + r[ f) a)+0(a 3 ) (j=l,2) 



S 



2 a 



(21) 
(22) 



Applying to this TPS of the sum a PA, say [1/1], and expanding the resummed result back in powers of the coupling 
a, we obtain 



5H/1] = 2 



2 a 



l- 5 (ri" 



(i) 



(2)n 

-r\ )a 



-(r\ 
2 y 1 



(1)2 (2)2 



2r 



(1U 2 ) 



)„ 2 



0(a 4 ) 



(23) 
(24) 



5 Formally, PA's are constructed for TPS's where the coefficients are independent of the expansion parameter. 
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The coefficient at a 3 of the power expansion of the result has a term 2r[ r\ , indicating that the result contains 
mixing effects at ^a 3 . This indicates some kind of interference effect between the two amplitudes for the observables 
(processes) S\ and S 2 - This is not acceptable because S is the sum of these two independent observables, i.e., the two 
contributions should be summed up incoherently. This argument remains basically unchanged when different PA's 
are applied, and/or when the TPS's are of higher order, either in a or in g s . 

It is true that expansion of Fe+m in powers of g s (using also expansion of toe in g s ) gives us a TPS of relatively 
high order (~ ~ a 5 ^ 2 ), and that we can thus apply PA's of relatively high order. Further, the higher order PA's 
are known to possess in general a weaker //-dependence than the corresponding TPS ]iS[ |. However, such an approach 
would predict nonphysical higher order effects and thus lead to unreliable predictions for the sum (Fe+m)- 

As argued in the previous sections, we follow a more conservative approach, by summing up by PA's each TPS (for 
Fe, we, Fm) separately. This approach gives us an additional freedom - to choose the RScl in each of these TPS's 
separately, in the natural range /x~ 1/r where r is a typical distance associated with each observable. For TPS's of 
Fg, toe, Fvi this implies the RScl choices: ue^oji (~2ttT); /i m ~mE; and /Jm^^e- Although, in this approach, 
we have to take as bases the TPS's (|l3|), (fl6|), ( fl9| ) of very low order (NLO), and some of these TPS's show very 
divergent behavior, we will now show that the results nonetheless show remarkably weak dependence on fi and Ae- 



B. Numerical results 



In our numerical analysis we used, unless otherwise stated, for the number of active (massless) quark flavors nj = 3, 
and for the QCD coupling constant the reference value a s (M%, MS) = 0.334. We used the MS scheme, and for the 
/3-function we used, for definiteness, the [2/3] PA. 6 

In Fig. [I] we present the results for the screening mass toe as a function of the corresponding RScl fi m , for the 
temperature choice T = 1 GeV. All approximations are based on the (NLO) TPS (|l6|) . We see that the diagonal PA 
[1/1] (a) significantly reduces the RScl-dependence, in comparison to the LO and NLO TPS's, in accordance with the 
arguments of Ref. JL9|. Furthermore, the effective charge (ECH) method §9], ||0| |l| of fixing the RScl in the NLO 
TPS gives us a value (fixed by definition) not far from the PA [1/1] values /We will choose the "physical" screening 

mass to be determined by the condition m^/ 1 '(/x m )=/x m (=171%), which gives us the value = 1.9 GeV for T = l 
GeV. 

In Figs. ||(a) and |](b), we present the results for Fe and Fm as functions of the respective RScl's /ie and fiM- 
The factorization scale Ae was chosen, in the spirit of the hierarchy relations ( |Io[ ) , to be the logarithmic mean of the 

typical scales /ie ~ 2irT and /im ~ toe: Ae = \J 2nTrn^\ where to^ 1 ' =1.9 GeV is the aforementioned "physical" 
[1/1] (a) screening mass. The PA [1/1] (a) for F E is based on the NLO TPS We sec that the PA [1/1] for 

Fe has drastically reduced the RScl-dependence, and its values are very close to those of the ECH prediction. 

On the other hand, the situation with Fm is less favorable, as seen from Fig. || (b). The two PA's [1/1] (g s ) and 
[0/2] (gs) were constructed from the TPS in powers of g s (not a) of Eq. Jl9|), taking for toe, appearing in Eqs. ( fl8| ) 
and (|19l), the PA [1/1] (a) with the RScl's for to-e and Fm taken equal: /i m = (J,m- In the case of the (NLO) TPS 
result for Fvi, we used for toe the NLO TPS result as well. While the two PA's give reduced RScl-dependence of Fm, 
their values differ drastically. In the case of the PA summation of Fe and toe the diagonal PA's [1/1] (a) are more 
physically motivated than the off-diagonal ones [0/2] (a), due to the RScl-independence of [l/l](a(/i)) in the large-/3o 
approximation |l9|| . That is why we did not employ the [0/2] (a) PA's for Fe and toe- On the other hand, for Fvi we 



have a TPS in g s (lj) = 2-k^J a(\x) of Eq. ([l9|), and not in a(/i), and therefore the choice [1/1] (g s ) for the right-hand 
side (RHS) of Eq. ( Jl9| ) is not physically better motivated than the choice [Q/2](g s ). This can be seen also from the 
comparably weak RScl-dependence of both PA approximants in Fig. |^(b). In order to choose between the two, we 
have to study, in addition, their variation under the variation of the factorization scale Ae- This will be discussed 
just below. 

In Figs. ||(a) and |s|(b ), we present the AE-dependence of the approximants. The RScl's were fixed, in accordance 
with the hierarchies (|10|), in the following way: [1e = 2itT, /L«M = Mm =m E> where toe has the aforementioned "physical" 
value: toe = TOe' = 1.9 GeV, for the case of the PA approximants and 1.4 GeV for the case of the TPS's. Figure [3(a) 



6 As argued in Refs. J2l| , this PA gives us a reasonable quasianalytic continuation of the MS /3-function for values of a s (/i) up to 
Cs(p-) = it X 0.32 fa 1.0, i.e., for values where the TPS /3 gives highly u n rel i able results. 

7 A generalization of the method of diagonal PA's has been developed ]20l El|, which gives complete RScl-independence | B0j , or RScl- and 
RSch-independence ]2]| . However, we have here TPS's available only at the low NLO order, at which the aforementioned approximants 
basically reduce to those of the ECH method. 
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shows that the choice [1/1] (g s ) for the TPS of Eq. ( |l9| ) in Fm leads to strong, and thus unphysical, As-dependence 
of Fe+m- Thus we have to discard the [1/1](<? S ) result of the TPS ((l^). The choice [0/2](g s ) for this part leads, 
on the other hand, to a result for Fe+m which is remarkably stable under variation of Ae in the entire interval 
m$ < A E < 2ttT. 

Therefore, we conclude that the choice [l/l](a) for Fe and for toe, and the choice [0/2](g s ) for the TPS of Eq. jis|), 
i.e., the curve Fe[1/1]+Fm[0/2] in Fig. ||(a), is the least unreliable among all the curves. 

The aforementioned properties of various approximants under changes of the RScl's and of Ae remain qualitatively 
the same when we change the value of the temperature. In Figs. ||(a) and f|(b), we present the results for the 
approximant (Fe [1 / 1]+Fa/ [0/2] )/ Fideal as a function of the temperature T and of the coupling g s (2nT) [= 2n \J a(2-KT)] . 
The middle curve (mid A) is for the canonical choice of the scales, /xe = 2-kT and /xm = Mm — m E where mg is, for each 
T, determined in the aforementioned way: m^^ 1 \fj, m ) = fi m . The factorization scale is taken to be Ae = ^/JIeJ^u- 
The lower curve (low A) is for the choice of the scales being all at the lower extreme: ^e = MM = Mm = m E = Ae, 
where Toe is determined in the aforementioned way. The upper curve (high A) is for the choice of the scales at the 
upper extreme: /xe = Mm = Mm — A-E — 2-kT . We notice that the curve of the low A choice follows well the curve for 
the canonical (mid A) choice in the entire depicted region of T (of g s ). The other choice of PA's ([1/1] (a) for Fe and 
toe; [1/1] (gs) for Fm of Eq. (|l^)) gives values well outside this range - see Fig. ||(a). The TPS's [for Fe, to-e, and 
for Fm of Eq. dTaVI do give us values in qualitative agreement with the lattice results for a specific choice of scales 
(mid A) [see FigTg(a)], but the TPS values change drastically when some of these scales, in particular /xm, change 
see Fig. |(b)] . In Figs. |(a) and 0(b), we included, for comparison, the values of the TPS up to 0(g# ) in powers of 
g s {2irT), i.e., the sum of the TPS's ©-© and ©. We note that the values for the TPS up to 0{g A s ) would be 
above 1. 

As all the curves up to this point have been given for the choice of three active massless quark flavors (rif = 3), we 
present in Figs. ||(a) and||(b), the results for rif ~0, 2, 3, 4, 6 (and for the mid A choice of the scales). We see that the 
curve for rif = 4 differs little from that for rif — 3 in most of the covered parameter space, while the curve for rif = 
differs significantly. However, in the temperature region of particular interest (T ~ 0.1-10 GeV), the choices rif = 3 
or 4 are expected to be more realistic. 



C. Comparison with other approaches 



We will now compare the r esults presented in Figs. ^ and with those of some other approaches. 



As mentioned in Sec. |III A| , Kastening [|22j and Hatsuda |23j] applied high order PA's to the full TPS for the sum 
R = (Fe + Fm)/ Fdoai, i.e., s TPS in powers of <7 s (m) including ~ g^. Kastening p3] showed that, in the case of value 
g s (T) ps 1.1 [a s (T) = 0.1, rif = 6], the application of PA's [2/2}(g s ) and [2/3] (gT) reduced the RScl dependence in 
comparison with the TPS results. At larger values g s (2nT) w 2 [a s (2TrT) — 1/3, rif = 3] there was no significant 
reduction. On the other hand, Hatsuda fl23fl applied a kind of modified PA's, by postulating that the expressions in 
the numerator and the denominator of the PA's have zero coefficients for the term ~ g s . The RScl dependence was 
then shown to be significantly reduced for many such modified PA's ([2/3], [3/2], [2/4], [4/2]; with rif =4). Further, 
his curves for the R obtained as a function of a s (2i:T) were qualitatively similar to our curves in Figs. |](b) and ||(b), 
with i? m i n = 0.97-0.98. However, the differences between his c urves were significant, and there was no clear principle 



to choose any specific one of them. Furthermore, in our Sec. [II A we stressed a more physical point of criticism of 
this approach. 

A reliable comparison with lattice results is hampered by the fact that (QCD-)lattice calculations have reproduced 
thermodynamic quantities (free energy or pressure, entropy, etc.) only for temperatures between the critical temper- 
ature T c and 4.5T C (i.e., for T < 1 GeV). In this temperature region the resulting values for R are less than 0.87 when 
rif — [Q, significantly lower than 1. If the number of (massless or light) quark flavors rif is larger than zero, the 
finite cutoff effects are not quite under control and are estimated (|] to increase the calculated R{rif) by about 15%, 
giving at T = 3.5T C (w 0.7 GeV) the values R w 0.90 ± 0.04 for n f = 2 (cf. Fig. 3 of Ref. ||). On the other hand, our 
approach, which is based on improved perturbative considerations, is expected to be reliable only for much higher 
temperatures T > 10 GeV where the corresponding values of the effective coupling parameters g s (27rT) are not much 
higher than 1. At such high temperatures we predict R w 0.975 or higher, for rif < 4 [cf. Fig. 0(a)], i.e., very near 
to 1. These predictions are not incompatible with the lattice results, however. For example, if we take from Fig. 3 of 
Ref. [S] the continuum estimate for p/T 4 (with rif = 2), which is given as a band-curve for T < 3.5T C (« 0.7 GeV), 



We note that in the thermodynamic limit F = —p, where p is pressure. Thus, R = p/psB where psB=Pidoal is the ideal gas pressure. 
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and extrapolate it smoothly to T > 10 GeV, we can obtain a fair agreement with our results. Further, we note that 
for T > 10 GeV, g s (2irT) can still be rather high \g s (2nT) = 1.1-1.2] and the numerical evaluations of perturbative 
TPS for such g s (2irT) differ from our results significantly. For example, the TPS up to 0(g 5 ), partly visible in Fig. ||, 
gives us for 3 s (2ttT) ps 1.2 (1.0) [o T « 10 (100) GeV] the value 0.925 (0.952). 

One interesting point in connection with lattice data is the dependence of R on the number nf of light or massless 
quark flavors. It is well known that the ideal gas expression for F shows an increase for |F| with increasing n/ , Eq. (|l^). 
Further, while lattice calculations show that the ratio R(rif) = F(n/)/F(n/)ideai also increases with increasing n/ at 
T ^T c , they indicate an inversion of the n/-dependence of R at the highest available T values (T w 4T C ) This is 
in accordance with our finding [see Fig. ||(a)] that at T ^> T c the ratio R decreases with increasing n/. 

A calculation in the <I>-derivable approximation scheme using hard thermal loop (HTL) propagators |l5| gives lower 
minimal values of R (i? m ; n ~ 0.82, see Fig. 3 there, where rif = 0), and the minimum is at a much higher value 
of g s (g s w 2.5) than our curves [g s {2 / nT) ss 1.1]. However, the <I>-derivable approximation was performed at the 
leading- loop order, and the expansion of the result R in powers of g s underestimates the positive term by a factor 
of 1/4. This indicates that the correction of this effect would push the results for R higher. 



IV. RESUMMATION RESULTS IN THE CASE OF </> 4 THEORY 



We can apply the same methods of resummation in the massless scalar 
corresponding TPS results from Ref. H. The interaction is 



Ant — 



4! v 



theory. In this case, we can use the 



(25) 



and the /3-function in the MS scheme is known to five loops |3 

da „ 9 17 



d In ji 



3a z 



32.54a 4 - 271.6a 5 + 2848.6a 6 + C(a 7 ) , 



(26) 



where a(/x) = g 2 (fi) /(16ir 2 ). For definiteness, we choose the PA [3/3] for this /3-function, in order to simulate better 
the running in the l&Tge-g(p.) region. The high energy contribution F$ to the free energy F is 



Fs(T) 



90 



T 1 



^s(T) 



^ 2 ){ 



1 



3 In 



4ttT 



9.29324 



i(» 2 )}+0(a 3 ) • 



(27) 
(28) 



In contrast with the analogous QCD quantity Fe, Fs does not depend on the factorization scale A^ (in QCD: Ae) at 
the available order. The TPS for the Debye screening mass ttil is 



m|(A F ) 



" 2tt 2 T 2 



ml (Af) 



-3 In 



4ttT 



\AttT 



5.39289 



(29) 



In contrast with the analogous QCD quantity ttje, the mass ttil does depend on the factorization scale. The low 
energy contribution Fl to F has the following associated TPS: 



12tt Fj 
T ml(A F ) 



= -1 



3 J_ 

2 4^ 



5(M 2 ) 



4 In 



A, 



4m L (A F ) 



3 1 



2 16tt 



r3 2 (M 2 ) 



0(g 3 ) . 



(30) 



The factorization scale dependence now cancels out within Fl, up to the available order (~ g 5 ). We have the hierarchy 
( |Io| ) as in QCD, with the substitutions itie t— * ttil and Ae i— > Ap. The entire analysis of the QCD case, as described 
in the previous section, can now be repeated following the same procedures. For the coupling parameter we chose 
the reference value g(2n GeV) = 4. For the temperature choice T = 1 GeV, the results for the screening mass 
as a function of the corresponding RScl fx m are presented in Fig. ||, and for Fs and FL as functions of the 
corresponding RScl's /is and in Figs. 0(a) and @(b). The factorization scale Af was determined similarly as in 

27rTm^°\ where in turn rnf^ is the "physical" screening mass determined by the condition 

: 0.882 GeV and A F = 2.354 GeV (for T = 1 GeV). 



the QCD case: A F ■ 
m£ (/i m ;Ajr) = n m (= 7Ti^ ) , which gives us the values m 



(0) 
L 
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The PA curves for in Fig. uYb) were obtained by applying the corresponding PA's [1/1] (g) or [0/2] (g) to the RHS 
of Eq. ( |30| ) and employing for m L the PA [l/l](a) at the same RScl (a m = ^ L ). The TPS curve was obtained by 
using the TPS of the RHS of Eq. @, using for m L the (NLO) TPS @. The [0/2] curve in Fig. 0(b) is much less 
RScl-dependent than [1/1], but not significantly less than the TPS curve. 

The dependence of the results on the factorization scale A^ is depicted in Fig. |[ We see that the PA choice [l/l](o) 
for Fs (and itil) and [0/2] (g) for Fl results in weak A^-dependence which is comparable with the TPS result, 9 while 
the [1/1] {g) choice again gives unacceptably strong A^-dependence. 

When the temperature is changing, so is g(2irT) according to the RGE (26) (we use [3/3] PA on the RHS). In 
Fig. ||(a) we present the results for the choice [1/1] (a) (for F s and mi) and [0/2](</) [for the RHS of Eq. (|o|)], and 
for the NLO TPS choice [for F§, mi, RHS of Eq. (|3(j)]. The "middle" values (mid A) are chosen for the relevant 
scales /is = 2-7rT, /LtL = /i TO = m£ 0u m ; Ap), and Ap = ^//^sMl- If these scales are changed to low A [fj, s =fx h =fj, m =Ap 

with m^ 1 ^ (Mm! Ap) =^t m ], or to high A (ns=HL = fi m = Ap = 2irT), the PA results vary with g(2irT) as presented in 
Fig. |(b). 



V. SUMMARY AND CONCLUSIONS 



We have presented here a method to improve the predictive power of perturbative thermal field theory. Predic- 
tions obtained by ordinary (truncated) perturbation theory (including a resummation to get rid of finite-T infrared 
divergences) suffer from serious divergence and renormalization scale (RScl) ambiguity problems. Therefore, a careful 
reorganization of these series is needed, which converts their physical content into expressions with better convergence 
behavior and greater stability under the variation of the RScl. 

Efforts in this direction have been undertaken by various authors during the last five years and considerable 
improvement has been achieved. However, most of these attempts have improved either the divergence problem or 
the RScl-ambiguity, but not both. 

We suggest an approach that addresses both problems, but concentrating mainly on the unphysical RScl- 
dependence. It is based on a physically motivated separation of the free energy density F into parts which are 
separately RScl-independent since each of them has an empirical significance on its own. This can be achieved 
because the separation principle is determined by observable effects (screening of the different massless degrees of 
freedom). In this way one simultaneously obtains a clear separation between terms which include (ring-diagram) 
resummation effects on the one hand, and purely perturbative contributions on the other hand, which are free of 
any such resummation. This is gratifying since it allows us to consistently apply Pade approximants (PA's) to the 
available truncated perturbation series (TPS's) of the individual parts, thereby avoiding the danger of double counting 
and spurious interferences of contributions from different kinematical regimes. We also Pade resum the (observable) 
screening mass which appears in the long-distance part of the free energy density. The resulting expression for the free 
energy density has strongly reduced RScl-dependence. In addition, we expect this method to show a good convergence 
when applied to TPS's of higher order, based on the well-known behavior of PA's p8| . 

Due to its close connection with physical effects (screening), we consider our approach to be less ad hoc and more 
physically motivated than some of the previous methods. 

As a consequence of the aforementioned separation of the free energy into (two) parts, the underlying TPS's for 
the construction of PA's were of low order. Regarded from a purely numerical point of view, this should alarm us and 
we should expect significant (nonphysical) instabilities of the resulting low order PA's under variation of the RScl fi, 
and instabilitcs for the sum of the two parts under variation of the factorization scale. However, our results both in 
QCD and in </> 4 theory are remarkably stable under both variations. This confirms additionally that the described 
separation forms a sound basis for the application of PA's. Further, thus resummed values of R = -F/-F]deai turn out 
to be below the value 1 for a relatively wide interval of the coupling parameter: g s (2irT) < 1.87 [<=> a s (2nT) < 0.278], 
i.e., T > 0.4 GeV, in QCD with n/ = 3; see Fig. f|(b). However, the method probably breaks down already at lower 
g s (2nT) » 1.0-1.2, i.e., T » 10-100 GeV, where it gives the local minimum i? min » 0.977 (if n f = 3). The results at 
T > 10 GeV do not contradict the results of lattice calculations. 



The latter is Ap'-dependent because we apply in Eq. (BOj) in this case the (NLO) TPS value for mj . 
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FIG. 1: The Debye screening mass as function of the rcnormalization scale (RScl) /i m , when T=l GeV. The upper of the two curves 
LO TPS and NLO TPS, respectively, has a(/i^) evolved by the one-loop and two-loop RGE from a(m^). All the other curves have a(/i^) 
evolved by the four-loop PA [2/3] beta function. 
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FIG. 4: The approximant .Re[1/1]+Rm[0/2] as function of temperature (a) and of g s (27rT) (b), for three different choices of the relevant 
scales (/ie, /im =fJ. m , Ae), as explained in the text. For comparison, the TPS curve in powers of g 3 (27rT) is included. In our case we have 
g 3 (27rT) « 1.2 (1.0) for T = 10 (100) GeV. 
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RScFs. Other details are given in the text. 
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FIG. 8: Various approximants for the normalized sum J?s+L = ^S+lV-Fideal as functions of the factorization scale Ap. Other data given 
in the figure and in the text. 




FIG. 9: The normalized sum i?s+L = -Fs+L/-Fideal as a function of g(2irT) (a), (b) presents variation of the rcsummcd result ijg[l/l] + 
-Rl[0/2] for three choices of the relevant scales (fi<j, f*L = /J-m, Ap), as explained in the text. 



